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We study the quantum effects of radiation pressure in a 
high-finesse cavity with a mirror coated on a mechanical res- 
onator. We show that the optomechanical coupling can be de- 
scribed by an effective susceptibility which takes into account 
every acoustic modes of the resonator and their coupling to 
the light. At low frequency this effective response is similar to 
a harmonic response with an effective mass smaller than the 
total mass of the mirror. For a plano-convex resonator the ef- 
fective mass is related to the light spot size and becomes very 
small for small optical waists, thus enhancing the quantum 
effects of optomechanical coupling. 

PACS : 42.50.Dv, 42.50.Lc, 43.40,+s 

Radiation pressure exerted by light plays an important 
role in quantum limits of very precise optical measure- 
ments. Quantum noise in interferometers has two funda- 
mental sources, the photon noise of the laser beam and 
the fluctuations of the mirror's position due to radiation 
pressure. Both lead to a quantum limit for measurement 
sensitivity and potential applications of squeezed states 
have motivated a large number of works in quantum op- 
tics @,|. 

When a movable mirror is exposed to a laser beam, 
its position is coupled to the laser intensity via radia- 
tion pressure ||. This optomechanical coupling may be 
enhanced using a high-finesse optical cavity that is very 
sensitive to small mirror displacements. Such a device is 
equivalent to cavities containing a Kerr medium and it 
has been studied for squeezing generation Q-|7j or quan- 
tum nondemolition measurements fa-jOl. 

Theoretical calculations of these effects are usually 
based on a description of the mechanical motion as a 
single harmonic oscillator. This model is well adapted to 
the motion of the center of mass for suspended mirrors, 
such as the ones of gravitational wave antennas p] , ^2[ . 
Radiation pressure can however excite internal acoustic 
modes of the mirror for which a simple description as 
a harmonic oscillator is not appropriate. These internal 
vibrations induce a deformation of the mirror which can 
be coupled to the light. 

In this paper, we study the effect of optomechanical 
coupling on quantum fluctuations when internal acoustic 
modes of the mirror are considered. We show that the 
coupling can be described by an effective response of the 



mirror which takes into account all acoustic modes and 
their spatial matching with the light. Similar results have 
been obtained for the Brownian motion of mirrors which 
have been studied for gravitational-wave interferometers 
|l3| p^| . We show that the same effective response de- 
scribes both thermal and radiation pressure effects. This 
response can be approximated to a harmonic response at 
low frequencies and the effective mass of this equivalent 
pendulum can be much smaller than the total mass of the 
mirror. This small effective mass enhances the quantum 
effects of optomechanical coupling. 

The system studied in this paper is sketched in figure ^. 
It consists of a single-port cavity with a movable mirror 
coated on the plane side of a mechanical resonator. The 
geometry of the resonator determines the spatial struc- 
ture of the acoustic modes. The model presented in this 
paper is however valid for any geometry. We first recall 
the basic properties of the optomechanical coupling us- 
ing a simple description in which both the field and the 
mechanical resonator are treated as one-dimensional ob- 
jects and the mirror motion is described as a harmonic 
oscillator (Section ||). 

We then study the effect of internal acoustic modes by 
taking into account the spatial structure of both the light 
and the resonator. We first determine the effect of a de- 
formation of the resonator on the light field (Section ^|). 
We then study the mechanical motion of the resonator 
when its plane side is submitted to the radiation pressure 
of the intracavity field. Wc finally define an effective sus- 
ceptibility which describes the optomechanical coupling 
between the light beam and the resonator (Section III). 
To illustrate the quantum effects of radiation pressure we 
study the quantum noise reduction of the field reflected 
by the cavity (Section [Tv| ). 
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FIG. 1. Model of a Fabry-Perot cavity with a movable mir- 
ror. The mirror is coated on the plane side of a plano-convex 
mechanical resonator 

We apply these results in the last section to the case 
of a plano-convex resonator (Section [v]) . Acoustic modes 
are then confined near the central axis of the resonator 
and their spatial structure can be described by analyti- 
cal expressions which are quite similar to gaussian optical 
modes of a Fabry-Perot cavity. We show that this geom- 
etry leads to a drastic reduction of the effective mass. 
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I. ONE-DIMENSIONAL MODEL OF 
OPTOMECHANICAL COUPLING 

We recall in this section the basic properties of the op- 
tomechanical coupling. For this purpose we neglect the 
internal motion of the resonator and we assume that the 
mirror moves without any deformation of its surface. The 
light is only sensitive to the mirror motion in the z di- 
rection (Fig. [l]) and this motion can be completely char- 
acterized by the position z (t) of the mirror at time t. 

The effect of the mirror motion on the intracavity field 
is a phase shift ip related to the change of the optical 
path followed by the light 



ip(t) = 2kz(t) 



(1) 



where k is the wavevector of the light. The cavity detun- 
ing thus depends on the cavity length and couples the 
field to the mirror motion. 

For small mirror displacements, the linear response 
theory shows that the Fourier transform z [f2] of the mir- 
ror position is proportional to the applied force Jl7[ 



z [Q] = x [Q] (2hkl [Q] + F T [Q]) . 



(2) 



where x is the mechanical susceptibility of the mirror. If 
we assume that the mechanical motion is harmonic, this 
susceptibility has a Lorentzian shape 
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(3) 



where M is the mass of the mirror and Qm its resonance 
frequency. The loss angle $ [f2] characterizes the damp- 
ing of the motion and is related to the quality factor Q of 
the resonance by $ [CIm] — l/Q- The first force in (||) rep- 
resents the radiation pressure exerted by the intracavity 
field. It is proportional to the momentum exchange 2hk 
during a photon reflection and to the light intensity / nor- 
malized as the number of photons reflected on the mirror 
by unit time. The second force Ft is a Langevin force 
describing the coupling of the resonator with a thermal 
bath. Its spectrum St [fJ] is related to the mechanical 
susceptibility by the fluctuation-dissipation theorem |17f] 



StM 



2k B T 



Im 
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xM 



(4) 



where T is the temperature and k B the Boltzmann con- 
stant. 

Equation (Q) shows that the radiation pressure exerted 
by the light couples the mirror motion to the light inten- 
sity. The mean effect of this coupling is a mean displace- 
ment z obtained from equation (^) by a statistical average 
of z [O = 0]. From equation (|l|) this leads to a non-linear 
phase shift for the mean field in the cavity equal to 



where I is the mean intracavity intensity. This intensity- 
dependent phase shift is equivalent to an optical Kerr 
effect. The cavity with a movable mirror is thus equiv- 
alent to cavities containing a Kerr medium which have 
been studied for squeezing generation 18-20] or QND 
measurements 21 2 2f| . The main differences with a pure 
Kerr medium are due to the dynamics of the moving 
mirror characterized by the frequency dependence of the 
susceptibility \ [^] an d to the presence of thermal noise. 

The non-linear phase shift i/jnl is an important pa- 
rameter to determine the efficiency of the optomechani- 
cal coupling. Quantum effects are significant if this phase 
shift is of the order of the cavity losses jl(J ■ In this case 
the displacement z induced by the mean radiation pres- 
sure is of the order of the width A/ T of the cavity res- 
onance (A is the optical wavelength and T is the cavity 
finesse). This condition actually corresponds to the ob- 
servation of bistability in the cavity. As usual in quantum 
optics, quantum effects are important near the bistabil- 
ity turning points. This condition depends on the optical 
characteristics of the system, such as the cavity finesse 
or the light intensity. It also depends on the mechan- 
ical properties of the resonator (eq. . In particular, 
quantum effects induced by optomechanical coupling are 
inversely proportional to the mass of the resonator. 



II. LIGHT REFLEXION ON A MOVING MIRROR 

We have shown in the previous section that the op- 
tomechanical coupling is based on two complementary 
effects. The first one is the phase shift of light induced 
by the mirror motion. The second effect is that the mir- 
ror moves in response to radiation pressure. We examine 
in this section the first effect when internal modes of the 
resonator are taken into account. 

Light is sensitive only to longitudinal displacements of 
the mirror. A longitudinal deformation of the plane side 
of the resonator can be described by its displacement 
u(r,z — 0, t) in the z-direction at every radial point r 
of the surface (the origin of the cylindrical coordinates 
is taken at the center of the mirror, see figure |l|). The 
electric field in front of the moving mirror can be written 
as 



E(r,t)=v (r)a(f)e""°', 



(6) 



where a (t) is the slowly varying amplitude of the field, loq 
the optical frequency and vq (r) the spatial structure of 
the beam in the z = plane. In the paraxial approxima- 
tion, the optical modes of the cavity are gaussian modes 
with their waist at position z = (23). Assuming the 
incident beam matched to the fundamental mode of the 
cavity, the spatial structure vq (r) of the field is given by 



vPnl = ±hk 2 x [0] J, 



(5) 



va (r) = 




(7) 
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where wq is the optical waist which depends on the ge- 
ometry of the cavity. 

The spatial structure of the reflected field is modified 
by the mirror motion since the field wavefront is dis- 
torted and reproduces after reflection the shape of the 
mirror. At every radial point r of the mirror, the opti- 
cal path followed by the light is changed by the displace- 
ment and the field experiences a local phase-shift equal to 
2ku (r, z = 0, t) . The reflected field E' (r, t) is then equal 
to 

E' (r, t) = v (r) a (t) e ^u(r,z=v,t) _ ^ 

Noting {v n (r) } the basis of gaussian modes of the cavity, 
one can write the reflected field as a sum over all modes 

E' (r, t) = (voe^^^vn) v n (r) a (t) e~^\ (9) 

n 

where the brackets stand for the overlap integral in the 
z = plane 

(/,<?>= / d 2 rf (r) g (r) . (10) 

Jz=0 

Equation (^) shows that the mirror deformation induces a 
diffusion of the light into all optical modes. This diffusion 
is however limited by the cavity and it becomes negligible 
for a non-degenerate and high-finesse cavity. In this case, 
modes in the sum ^ evolve at a frequency about the res- 
onance frequency luq of the fundamental mode whereas 
the differences between the resonance frequencies of these 
modes are large compared to the cavity bandwidth. As a 
consequence, all modes except the fundamental one are 
filtered by the cavity bandwidth and cannot propagate 
in the cavity. One can show that the diffusion in those 
modes is equivalent to losses for the fundamental mode 
and that these losses become negligible for a high-finesse 
cavity | p4| . Only the fundamental mode has thus a sig- 
nificant contribution in the sum (|^) and we obtain 

E' (r, t) = {y e 7ik ^' =0 ' t \vo) E (r, t) . (11) 

For small displacements u this expression can be approx- 
imated to 

E' (r, t) w [1 + 2ik (u (z = 0, t) , v 2 )] E (r, t) . (12) 

The field thus experiences a global phase-shift tp (t) which 
can be written as 

ip{t) = 2ku{t), (13) 

where u (t) is the displacement of the mirror averaged 
over the optical waist 

u(t) = (u(z = 0,t),v 2 ). (14) 

The effect of the mirror motion on the intracavity field 
is thus equivalent to the one obtained in the one- 
dimensional model. It corresponds to a phase shift given 



by an equation similar to (|lj) where the one-dimensional 
displacement z (t) is replaced by the averaged displace- 
ment u(t). All results concerning the effect of the mirror 
motion on the field can thus be generalized to a mechan- 
ical resonator. The light is only sensitive to the displace- 
ment u of the mirror which takes into account the spatial 
overlap between the intracavity field and the mirror mo- 
tion. In the next section we determine this displacement 
when the resonator is submitted to the radiation pressure 
of the intracavity field. 

III. RADIATION PRESSURE EFFECTS 

Spatial and frequency characteristics of the acoustic 
modes depend on the geometry and acoustic properties 
of the mechanical resonator. To determine the effect of 
radiation pressure on the mirror motion, it is only nec- 
essary to assume that these modes are described by a 
set of displacements {Tt n i^)} which forms a basis of 
the resonator motion. Each mode n obeys a propagation 
equation inside the resonator [p5[ 

c 2 t ATt n + (cf - c\) ^ (^.Tt n ) + fi£"t?„ = 0, (15) 

and fulfills the boundary conditions corresponding to a 
free resonator (see Appendix) . In this equation Q n is the 
eigenfrequency of mode n and cj, Ct are the longitudinal 
and transverse sound velocities related to the Lame con- 
stants A and /i of the resonator and to its density p by 

ci = ^(X + 2p) /p, (16a) 
c t = \J pip- (16b) 

Any displacement Tt (T*, t) can be expressed as a linear 
combination of the acoustic modes T? n (~f*) 

it t) = a„ (t) lt n (J>) , (17) 

n 

where a n (t) is the time-dependent amplitude of mode n. 
Using this decomposition we can determine the evolution 
of it t) when a radiation pressure force is applied on 
the plane side of the resonator. The total energy E of 
the resonator is the sum of the kinetic energy, the poten- 
tial energy and the energy associated with the external 
force. It can be decomposed in the following form (see 
Appendix) 




M n represents the mass of the acoustic mode n which is 
proportional to the volume of the mode inside the res- 
onator 
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M n =p I d s r 
lv 



(19) 



The term (l? 



id j " 



in equation ( |18[ ) represents the 
spatial overlap of the scalar product between the radia- 
tion pressure ~P r ad and the acoustic mode Tt n (eq. 10). 
The radiation pressure is directed along the z-axis and 
its amplitude F ra< i (r, t) at radial position r of the plane 
side and at time t is related to the intracavity intensity 
I (t) = \a(t)\ 2 and to the spatial structure vq of the field 
by 



F rad (r,t)=2hkl(t)vl (r) 



(20) 



The total energy (eq. |l8|) appears as the sum over all 
modes n of the energies of forced harmonic oscillators. 
From Hamilton's equations one deduces that each mode 
amplitude a n (t) obeys the evolution equation 



d a Tl 



+ Sl 2 n a n 



idi It ', 



dt 2 ' M, 
This equation can be written in the Fourier space as 



(21) 



a n [St] = Xn [St] {Prad [SI] 



(22) 



where Xn [St] = 1/M n — tt 2 ) is the susceptibility in 
absence of dissipation for a harmonic oscillator of mass 
M n and eigenfrequency fl n . The resonator can thus be 
considered as a set of independent harmonic oscillators, 
each oscillator being associated with an acoustic mode. 
These oscillators are driven by an external force which 
corresponds to the projection of the radiation pressure 
onto the spatial structure of the acoustic mode in the 
z = plane. 

Up to now we have assumed that the resonator has 
no damping. The coupling with a thermal bath can be 
deduced from Navier-Stokes equation [p5| or from a gen- 
eralization of the approach used in the one-dimensional 
model. Each acoustic mode is indeed equivalent to a har- 
monic oscillator and its damping can be described by a 
dissipative part added to the mechanical susceptibility 
and by an additional Langevin force. The susceptibility 
Xn of the acoustic mode n has thus an expression similar 
to the one-dimensional case 



Xn [St] 



1 



M n (02 - W _ zfi 2 $n [n] y 



(23) 



where <&„ [tt] is the loss angle of mode n. The amplitude 
a n of mode n is now given by 



[fi] = Xn M ((frad M 



F T ,n M 



(24) 



where Fr, n is a Langevin force describing the coupling 
of mode n with the thermal bath. We assume that these 
forces are statistically independent from each other and 
that their spectra are related to the susceptibilities Xn 



through fluctuations-dissipation theorem (eq. [|). Acous- 
tic modes are then independent and the equation of mo- 
tion for each amplitude a n (eq. |24|) corresponds to the 
usual expression for a damped harmonic oscillator driven 
by the projection of the radiation pressure. 

We have shown in the previous section that the effect 
of the mirror deformation on the intracavity field only 
depends on the longitudinal displacement u (t) which cor- 
responds to the displacement in the z = plane averaged 
over the beam waist (eq. |lj). From equations (|l4|), (|lj), 
( pp| ) and (^4|) this displacement can be expressed in terms 
of the intracavity intensity / and of thermal fluctuations 



:[fi\ = X eff p] (2hki [ci] + F T [ii]) 



(25) 



Xef f appears in this equation as an effective susceptibility 
given by 



X 



(26) 



where u n stands for the z-component of the acoustic 
mode ~it n . The effective susceptibility Xeff is then equal 
to the sum of all susceptibilities Xn weighted by the over- 
lap between the acoustic mode and the transverse inten- 
sity distribution Vq. The force Ft in equation (|25|) is an 
effective Langevin force related to the forces Fr, n of each 
acoustic modes by 



F T [SI] 



Xn M 
Xeff [SI] 



F T> n [SI] 



(27) 



The spectrum of Ft can be determined using the inde- 
pendence of the Langevin forces Ft. u and their relation 
to susceptibilities Xn ( ec l- %>■ We find that the force F T 
is related to Xef f by the fluctuation-dissipation theorem 



S T P] = ^— Im 



1 



Xeff [SI] 



(28) 



This means that the resonator in absence of external 
force is in thermodynamic equilibrium at temperature T. 
Equation also shows that the effective susceptibility 
Xef f describes both the effects of radiation pressure and 
of thermal noise represented by the Langevin force Ft- 

Results obtained here are similar to the ones presented 
in section || for a one-dimensional model. The longitudi- 
nal displacement u is related to the intracavity intensity 
/ and to thermal fluctuations by an expression similar to 
equation (^|). The coupling between the resonator and 
the gaussian light beam is completely described by the 
effective susceptibility Xeff- This susceptibility indeed 
takes into account all acoustic modes of the resonator and 
their spatial matching with the light. Treatments made 
in the framework of the one-dimensional model can thus 
be generalized by replacing the mirror displacement z (t) 
by the longitudinal displacement u (t) and the harmonic 
susceptibility x by the effective susceptibility Xef f ■ This 
is illustrated in the next section where we study the quan- 
tum noise reduction of the field reflected by the cavity. 
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IV. QUANTUM-NOISE REDUCTION 

Squeezed-state generation has already been studied for 
a single-ended cavity containing a pure Kerr medium 
pTsj pti[ or for a cavity with a harmonically suspended 
mirror S. In this section we extend these results to the 
case of a mirror coated on a mechanical resonator. We 
thus consider the system sketched in figure [l] and we de- 
termine the quantum fluctuations of any quadrature of 
the reflected beam. 

For a nearly resonant high-finesse cavity, we have the 
following relations between the complex amplitudes a m , 
a and a out of the incident, intracavity and reflected fields 
respectively 



[7 - »* (t)] a (t) + v/2 7 a m (t) 



da 
~dt 



(29a) 
(29b) 



The first relation determines the dynamics of the intra- 
cavity field a. t is the round trip time, 7 the damping 
rate of the cavity (1 — 7 and \f2Fj are respectively the re- 
flection and transmission of the input mirror, with 7 <C 1) 
and ^ is the cavity detuning assumed small compared to 
1. $ depends on the cavity lentgh and couples the field 
to the mirror motion. It can be written as the sum of the 
cavity detuning Wo without light and the effect of the 
averaged displacement u (t) (eq. Oh 



* = *o + 2Kt. 



(30) 



The mean cavity detuning W differs from Wo by the 
contribution due to the mean mirror displacement. Using 
equation (|2^) this contribution appears as a non-linear 
phase shift W jv l related to the mean intracavity intensity 
I by an expression similar to equation (^|) 



y NL = 4hk 2 Xeff [0]T. 



(31a) 
(31b) 



The behavior of the mean fields is thus the same as the 
one for a cavity containing a pure Kerr medium. In par- 
ticular the non-linear phase shift is responsible for the 
bistability of the system ||. From equations ( p9| ) the 

mean intracavity intensity is related to the incident power 
pin by 



( 7 2 + I = 2 7 Ap»\ 



he 



(32) 



Due to the intensity dependence of the non-linear phase 
shift, one can obtain different intracavity intensities for 
a given incident power. This bistable behavior can be 
characterized by the inverse slope a = dP m /dI of the 
curve giving the intracavity intensity as a function of the 
incident power, a is proportional to 



The positive, negative, and null values of a correspond, 
respectively, to stable branches, unstable branch, and 
turning points pp| ]. 

To determine the fluctuations Sa out [ft] of the field re- 
flected by the cavity, we use the semiclassical method in 
which quantum fluctuations are treated as classical ran- 
dom variables associated with the Wigner distribution 
pp| p8| . Their evolution is deduced from the classical 
equations ( |29| ) linearized around the mean state. We ob- 
tain 



(7 - iW - iftr) 5a [ft] = ^/2^5a m [ft] + ia5^ [ft] , (34a) 
5a out [ft] = y^Sa [ft] - Sa ln [ft] , (34b) 
atfW [ft] = Xeff [ft] ^nl (6a [ft] + 5a* [ft]) 
+2ka Xe ff [ft] F T [ft] , (34c) 

where a is the mean intracavity field and Xef f [ft] = 
Xeff [ft] I Xeff [0] is the mechanical susceptibility normal- 
ized to 1 at zero frequency. We deduce from these equa- 
tions the input-output relations for the field which give 
the output field fluctuations 5a° ut as a function of the 
input field fluctuations 5a m and of the Langevin force 
Ft 

5a out [ft] = {ci [ft] 5a m [ft] + c 2 [ft] 5a m * [ft] 

+c T [ft]F T [ft]}, (35) 

where the coefficients c\ [ft], c 2 [ft] and ct [ft] depend on 
the system parameters 

1 



d [ft] = - {(7 + iW) (7 + ^ + 2M NLX eff [ft]) 



(ftr) 2 } 



2% _ 

C 2 [ft] = -^l^NLXeff [ft] , 
21 



(36a) 
(36b) 



<7 oc 7 + * + 2^^ 



NL- 



(33) 



c T [ft] = ^ ^ka (7 + i* - iftr) X ef f [ft] , (36c) 

A = (7 - iftr) 2 + * 2 + 2^ NL Xeff [ft] • (36d) 

From these relations one can determine the spectrum 
gout f or an y q uac l ra ture ag ut of the reflected field defined 
by 

Sa° e ut [ft] = e- l6 5a out [ft] + e l6 5a out * [ft] , (37a) 
6a° e ut [ft] 5a° e ut [ft'] = 2tt5 (ft + ft') S° e ut [ft] , (37b) 

where the bar stands for the average over the Wigner dis- 
tribution. Equations ( ^B|) and ( ^6[) show that the quan- 
tum properties of the reflected field only depend on a few 
parameters. The cavity is characterized by the damping 
7 which is related to its finesse T = n/j, by the cav- 
ity bandwidth Q cav = 7/V and by the mean detuning ty. 
The mirror motion is described by the non-linear phase 
shift ^Pjv l and by the frequency dependence Xef f of the 
mechanical response. The last coefficient ct describes 
thermal effects associated with the Brownian motion of 
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the mirror. To obtain quantitative values for the noise 
reduction, it is necessary to determine the mechanical 
response of the resonator which may depend on its ge- 
ometry. We study in the next section the case of a plano- 
convex resonator. 

Note finally that the differences with a pure Kerr ef- 
fect are only due to the frequency dependence of the 
mechanical response and to thermal fluctuations. One 
gets the usual expressions for a Kerr medium by tak- 
ing Xef f [St] = 1 and Ft [Si] = in equations ( |35| ) and 



V. OPTOMECHANICAL COUPLING WITH A 
PLANO-CONVEX RESONATOR 

We consider in this section that the resonator has a 
plano-convex geometry with a mirror coated on its plane 
side (Fig. |l|). We first determine the analytical expres- 
sion of the effective susceptibility. We then study the 
quantum- noise reduction of the field reflected by the cav- 
ity and we finally derive the effective mass associated 
with the optomechanical coupling. 



A. Effective susceptibility 

If the resonator thickness ho is much smaller than the 
curvature radius R of the convex side, the propagation 
equation can be solved using a paraxial approxima- 
tion and one gets analytical expressions for the acous- 
tic modes corresponding to gaussian modes p^ j. Only 
modes that have a non-zero overlap with the light inten- 
sity contribute to the effective susceptibility (eq. ^6j). As 
a consequence we disregard shear modes which induce no 
longitudinal displacement and we consider only the com- 
pression modes that have a cylindrical symmetry. Those 
modes are defined by two integers, a longitudinal index n 
and a transverse index p, and the longitudinal displace- 
ment u n>v {r,z) at point of radial coordinate r and axial 
coordinate z is given by 



,(r,z) =e- r2/w *L p (2r 2 /w 2 n ) cos ( 



nw 
Mr)' 



(38) 



U n ,p is composed of a transverse gaussian structure with 
a waist w n , a transverse Laguerre polynomial L p and a 
cosine in the propagation direction, h (r) is the resonator 
thickness at radial position r given by 

Mr)«Ao-^- (39) 
Acoustic waists w n depend on the longitudinal index n 

^ = ^1^. (40) 

717T 



Acoustic modes solution of the propagation 

equation inside the resonator and each mode evolves with 
an eigenfrequency St ntP given by 



Si 



o 2 



i 2 + -\^n(2p+l) 

7T V H 



(41) 



where SIm = nci/ho- 

One can now derive an analytical expression for the 
effective susceptibility as an infinite sum over all modes 
{n,p} (eq. |26|). In this expression the mass M n of the 
acoustic mode {n,p} (eq. |l9|) only depends on the longi- 
tudinal index n and is equal to 



M„ 



-ph w ri 



(42) 



This mass is proportional to the volume of the acoustic 
mode and is smaller than the total mass of the resonator. 
We can also derive the overlap between acoustic and op- 
tical modes 



2wl 



2wi 



2wl + !ii„ V 2w 



(43) 



The overlap only depends on the ratio w n /wo between 
acoustic and optical waists. 

We assume for simplicity that the loss angles $ n ,p [St] 
are the same for all modes and are constant in frequency. 
As a consequence these loss angles are simply related to 
the quality factor Q of the fundamental mode {1,0} of 
the resonator by 



$„, p [St] = l/Q. 



(44) 



These equations allow to compute any parameters of 
the optomechanical coupling. For example th e non -linear 
phase shift ^ nl is related through equation ( 31b ) to the 
effective susceptibility at zero frequency which is given 
by 



Xeff 



n , p 



M n sil Z 



(45) 



This sum can be numerically computed. In the following 
we consider a small plano-convex resonator of thickness 
ho equal to 1.5 mm and with a curvature radius R of 
the convex side equal to 150 mm. For these values the 
acoustic waist w\ of the fundamental mode is equal to 3.8 
mm. It is much larger than the optical waist wo which 
is typically of the order of 100 /im for a cavity of 1 mm 
long with a curvature radius of the input mirror equal to 
1 m. As a consequence the overlap ( (43[ ) slowly decreases 
with the mode indexes {n,p\ and it is necessary to sum 
over more than 10 6 modes to obtain a precision better 
than 10 -3 . 
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B. Optimum squeezing 

We determine now the optimum squeezing S op t [ft] of 
the reflected field, that is the minimum value of the spec- 
trum S° g ut [ft] obtained at every frequency ft by scanning 
the quadrature 9. We assume that the mechanical res- 
onator is made of silica (density p = 2200 kg/m 3 , longi- 
tudinal sound velocity q = 5960 m/s). Its fundamental 
resonance frequency Qm/^ is then equal to 2 MHz and 
we take a quality factor Q of 10 6 . The cavity is char- 
acterized by its damping 7 equal to 10 -5 (cavity finesse 
T = 3 x 10 5 ) and its working point is defined by the mean 
detuning 'F chosen equal to — O.27. We also assume that 
the cavity bandwidth ft cav — 7/r is equal to Qm and 
that the optical wavelength A is equal to 800 nm. 

To determine the non- linear phase shift <F jvl 1 we com- 
pute the effective susceptibility Xef f [0] at zero frequency 
from equation ([45]). We obtain Xeff [0] = 1.4 x 10~ 8 m/N 
for an optical waist wq of 200 /im. From equations ( |31b| ) 
and (|3^) one then gets a non- linear phase shift ^nl equal 
to O.287 for an incident power P m of 10 mW. Note that 
these values correspond to a positive slope a of the bista- 
bility curve, equal to 0.937 2 (eq. |33|). The working point 
is thus on a stable branch of the bistability curve. 



sjai 



















\ 







Frequency Q. 

FIG. 2. Optimum noise spectrum Sopt [ft] of the reflected 
field as a function of frequency. Curves (a) and (b) are 
obtained at T = 4K and at zero temperature respectively. 
Dashed curve (c) corresponds to the pure Kerr effect for the 
same parameters ($ = — O.27, ^nl = O.287) 

Figure || shows the optimum squeezing S op t [ft] as a 
function of frequency. Curve (a) is obtained at a tem- 
perature of 4 K. The noise spectrum exhibits a strong 
reduction over a wide frequency range from zero up to 
the first resonance frequency of the resonator. For higher 
frequencies the spectrum shows an excess noise at every 
mechanical resonance frequencies. 

It is instructive to compare this spectrum to the ones 
obtained at zero temperature (curve b) and for a pure 
Kerr effect with the same parameters ^nl and 'F (curve 
c). At low frequencies (ft < Qm) the noise spectrum at 
zero temperature (curve b) is similar to the one obtained 
with a Kerr effect (curve c). The whole effect of the me- 



chanical motion can be interpreted as a non-linear phase 
shift for the light. At finite temperatures (curve 

a) , the Brownian motion of the mirror slightly increases 
the noise. Note that this effect depends on the dissipa- 
tion mechanisms in the resonator. In particular the large 
increase at zero frequency is due to the choice of a con- 
stant loss angle. One would obtain less thermal noise in 
the framework of a Navier-Stockes model for which the 
loss angle is a linear function of frequency. 

For higher frequencies (ft > Om) the dynamics of the 
mechanical resonator plays an important role. At zero 
temperature (curve b) we observe a series of dispersion 
shaped resonances centered on every mechanical reso- 
nance fi n ,p, with an important noise reduction for fre- 
quencies slightly below each resonance. Thermal noise 
however masks this behavior since the Brownian motion 
is concentrated around the mechanical resonances and 
increases the noise (curve a). 

To summarize this discussion, it appears that the most 
interesting frequency domain for quantum noise reduc- 
tion is the low frequency domain (ft < flj^)- For suffi- 
ciently low temperature the noise behavior is similar to 
the one obtained with a pure Kerr effect and it mainly de- 
pends on the non-linear phase shift >F jvi • More precisely 
the frequency dependence of the noise spectrum roughly 
corresponds to a low-pass filtering due to the cavity band- 
width and the amplitude of noise reduction becomes im- 
portant near the bistability turning points |2(J . To reach 
these points the non-linear phase shift 'Fat l must be large 
enough, that is of the order of 7. 

C. Effective mass 

The non-linear phase shift depends on the mechanical 
and optical properties of the system. In particular it is re- 
lated to the spatial overlap between the various acoustic 
modes and the light (eq. |I|). Figure || shows the optimum 
squeezing S opt [ft] for different values of the optical waist 
wq (400, 200 and 100 /im), the other parameters being 
identical to the ones of figure || (P,„ = 10 mW, 7 = 10~ 5 , 
* = -O.27, T = 4 K and Q = 10 6 ). For each curve we 
have computed the effective susceptibility Xef f [0] at zero 
frequency and we have determined the non-linear phase 
shift ^nl- Curves (a) to (c) correspond to increasing 
non-linear phase shifts (O.II7, O.287 and O.627, respec- 
tively) and to decreasing slopes of the bistability curve 
(0.997 2 , 0.937 2 and 0.797 2 , respectively). 

One observes a drastic increase of the noise reduction 
when wq decreases. This result can be interpreted as a re- 
duction of the effective mass associated with the optome- 
chanical coupling. Quantum effects observed in figures ^ 
and ^ in the low frequency domain are actually similar 
to the ones obtained with a harmonically suspended mir- 
ror. The effective susceptibility of the resonator is thus 
equivalent at low frequency to the susceptibility of a sin- 
gle harmonic oscillator with a resonance frequency flM- 
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The mass of this equivalent pendulum is however differ- 
ent from the total mass of the resonator since it depends 
on the spatial matching with the light. This effective mass 
M e f f can be deduced from a comparison between the ef- 
fective susceptibility Xeff [0] at zero frequency and the 
one of a harmonic oscillator (eqs. || and [f5|) 



1 



M P fftt 



eff lL M 



( u i u n,p) 

n.p 



v 2 



M fi 2 



(46) 



The effective mass is related to the masses of all acous- 
tic modes coupled to the light. The number of acoustic 
modes contributing to the sum strongly depends on the 
light waist. The typical variation length of the Laguerre 
polynomial L p (2r 2 /w 2 t ) in the spatial structure of the 
acoustic mode {n,p} is w n /y/p. The overlap (vo,u n ,p) 
is thus equal to 1 as long as wq is much smaller than 
w n/y/p (eq. [f3| ). For a given value of n, the number of 
transverse modes which contribute to the effective mass 
is proportional to wf/nw^ and it increases when wq de- 
creases. As a consequence the mass becomes smaller and 
quantum effects of radiation pressure become larger. 












a, 



Frequency £2 

FIG. 3. Optimum noise spectra S op t [fi] of the reflected 
field for different optical waists Wo. Curves (a) to (c) cor- 
respond respectively to 400, 200 and 100 /mi 



From equations (41) to (43) one can express the effec- 
tive mass as a function of the waists w\ , wq and of the 
mass Mi of the fundamental mode 



Mi 



4wf 



M 



eft 



n j p 



{2w\ + 



2w\ 



2w{ 



nuin 



2p 



n + 



k(2p+l) 



(47) 



Solid line in figure ^ shows the variation of the effective 
mass as a function of the ratio W\/wq for a resonator of 
thickness ho = 1.5 mm and of curvature radius R = 150 
mm. It clearly appears that this mass decreases with the 
optical waist and very small values can be reached. The 
mass Mi of the fundamental mode is equal to 37 mg (eq. 
E2I) and one gets an effective mass of 1 mg for an optical 



waist equal to 101/IO (380 /mi) and a mass of 0.2 mg 
for a waist of 100 /mi. Such small values would be very 
difficult to obtain with a harmonically suspended mirror 
for which the mass associated with the global motion is 
the total mass of the mirror. 
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FIG. 4. Effective mass M e ff of the resonator as a function 
of the ratio between fundamental acoustic and optical waists. 
Dashed curve represents the optical mass M opt related to the 
volume of the resonator illuminated by the light. Vertical 
scale is normalized to the mass Mi of the fundamental acous- 
tic mode (37 mg) 

One can get a simple physical insight into the effective 
mass in the case of a resonator thickness much smaller 
than the curvature radius of the convex side (ho <C R). 
We can then assume that the transverse modes of the res- 
onator are degenerate, that is the resonance frequencies 
£lnj> are independent of p. The sum over p in equation 
( [47] ) is a simple geometric sum and one gets an estimate 
M opt of the effective mass given by 



Mi 
M 



opt 



W 

2m 



2 1 

-Y — 



7T 

12 



m 

Wo 



(48) 



Using the expression of Mi (eq. [l2]) we finally obtain 



M, 



opt 



12 /7T 2 

^2 (jphowo 



(49) 



The term in brackets corresponds to the mass of the part 
of the resonator illuminated by the light beam. This op- 
tical mass is a good approximation of the effective mass 
M e f f as shown by the dashed curve in figure ^. We have 
thus shown that the effect of the resonator motion on 
the light is equivalent to the one of a harmonically sus- 
pended mirror of mechanical resonance frequency Vt^ 
and of mass related to the light spot size. This mass can 
of course become very small for a small optical waist. 



VI. CONCLUSION 

We have studied the quantum effects due to radiation 
pressure in a high-finesse cavity with a mirror coated on 
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a mechanical resonator. We have shown that the optome- 
chanical coupling between the gaussian laser beam and 
the acoustic modes of the resonator leads to a non-linear 
phase shift for the light. This phase shift is related to 
the intracavity intensity through an effective susceptibil- 
ity which takes into account all the acoustic modes and 
their coupling to the light. This susceptibility also de- 
scribes the effect on the light of the Brownian motion of 
the mirror. 

We have studied the quantum noise reduction of the 
field reflected by the cavity. This quantum effect mainly 
depends on the behavior of the non-linear phase shift at 
low frequency. In this frequency domain the mechanical 
response of the resonator can be approximated to a har- 
monic response. The effect of optomechanical coupling 
is then equivalent to the one obtained with a harmoni- 
cally suspended mirror of resonance frequency equal to 
the fundamental resonance frequency of the resonator. 
The mass of this equivalent pendulum is however smaller 
than the total mass of the mirror, thus enhancing the op- 
tomechanical coupling. We have shown that for a plano- 
convex resonator this effective mass is of the order of the 
optical mass which corresponds to the volume of the res- 
onator illuminated by the light. The effective mass can 
be two or three orders of magnitude smaller than the 
total mass of the mirror and large quantum effects are 
obtained for a reasonable input power when the optical 
waist is small enough. 

This drastic decrease of the effective mass seems to be 
a specific behavior of the plano-convex geometry. The ef- 
fective susceptibility of a cylindrical mirror has already 
been determined at low frequency for the study of ther- 
mal effects in gravitational-wave detectors 0Q. The 
reduction of the effective mass does not exceed a factor 
10 below the total mass of the mirror. The plano-convex 
geometry thus allows to get simultaneously a high res- 
onance frequency and a very small mass so that large 
quantum noise reduction can be obtained over a wide 
frequency range. One can of course obtain significant 
non-linear phase shifts with a cylindrical mirror but the 
fundamental resonance frequency is very low. 

Finally, we have only studied in this paper the quan- 
tum noise reduction of the field reflected by the cavity. 
Similar results would however be obtained for other quan- 
tum effects of radiation pressure such as the possibility 
of quantum nondcmolition measurement of light inten- 
sity jl(| or the quantum limit in interferometric measure- 
ments. For example it has been shown that the quantum 
limit for measurement sensitivity is proportional to the 
susceptibility characterizing the mirror motion |J. The 
quantum limit induced by internal acoustic modes is thus 
proportional to the effective susceptibility studied in this 
paper. In particular the sensitivity is reduced if mirrors 
have small effective masses. 
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APPENDIX A: PROPAGATION EQUATION FOR 
A RESONATOR SUBMITTED TO RADIATION 
PRESSURE 

In this appendix we derive the evolution equation of 
the acoustic modes for a resonator submitted to a radi- 
ation pressure force. This derivation is actually similar 
to the phonon decomposition in presence of an external 
force. 

The evolution equation of any displacement it 
can be deduced from the Lagrangian L = T — U where 
the kinetic energy T and the potential energy U are equal 

to 63 ' 



T 



1 



dt 



cfr, 



(Ala) 
(Alb) 



where p is the density of the resonator. The strain tensor 
u%j is given by 



1 (d 



2 V dxi 



dui 



(A2) 



where u$ (i = 1, 2, 3) are the cartesian components of 
the displacement it (T*,t) at point It of cartesian coor- 



dinates 



The stress tensor cr,,- is related to the strain 



tensor u^j by Hooke's law 



2fiUij + \ V.lt5i 



(A3) 



where A and /i are the Lame constants of the resonator. 
We deduce from the Lagrange equations the propagation 
equation of it (T*, t) in the resonator 



(A4) 



For a free resonator the constraint tensor cry at every 
point of the surface must satisfy the following condi- 
tion for any cartesian index i 



E 

3 



Oij t) rij = 0, 



(A5) 



where it is the normal vector at point T\ 

Acoustic modes are defined as monochromatic solu- 



tions lt n ("f > ) e~ lS1 "* of equations (A4) and (A5). Each 
acoustic mode can be decomposed into longitudinal and 
transverse components noted tt l n and it^ respectively 
and defined by 







0, ^. 



0. 



(A6) 



These components obey the following propagation equa- 
tions 



Alt 



o 2 

A U n + — U n 

c t 



(A7a) 
(A7b) 



where q and Ct are the longitudinal and transverse sound 
velocities (eq. [l6| ). 

Since the set {lt n (f^)} of aco ustic mode s fo rms a ba- 
sis for the solutions of equations ( |A4j ) and (|A5|), any dis- 
placement (~f*,i) can be decomposed on these modes 
with time-dependent amplitudes a n (t) (eq. ^) . We now 
determine the total energy E = T + U associated with 



this displacement. From equation (Ala) the kinetic en- 
ergy T is equal to 



T : 



^2 



M„ 



da„ 
~dt 



(A8) 



where M n is the mass of mode u n (eq . [19| ) . The potential 
energy U is given by equation ( Albj ) and can be written 
as 



U 



d((iijUj) dui 



dxi 



(A9) 



The first term can be transformed into a surface integral 
which is equal to zero for a displacement satisfying the 
boundary condition (A5). Using the decomposition of 
lt n into longitudinal and transverse components and the 
propagation equations (A7) one then gets 



U = E \ M ^i [o» (*)] s 



The total energy E is finally equal to 

2 



E 




1 



M n Q 2 n [a n (t)Y 



(A10) 



(All) 



and appears as the sum of the energies of free harmonic 
oscillators. 

Up to now we have assumed that the resonator is free 
from any external constraint. In presence of radiation 
pressure there is an additional contribution to the energy 
corresponding to the work of internal constraints opposed 
to the external force 



W = 



d 2 rPrad ("r*, t) .it (T\ t) . 



(A12) 



The total energy E is then equal to 

2 



E 



1 , , / da n 
2 Mn [-& 



\M n n 2 n [an (i)f 



^rad, Itn^ «„ (*)| , 



(A13) 



where (^Pradj^n) 1S the overlap integral in the z = 
plane of the scalar product between the radiation pres- 
sure and the acoustic mode u n (eq. 0). We have thus 
shown that the total energy is the sum over all modes 
of the energies of forced harmonic oscillators. The evo- 
lution equation ( pl| ) for each mode amplitude a n is then 
deduced from Hamilton's equations. 
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